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Abstract 

We study in detail the influence of molecular interactions on the Frank elastic con- 
stants of uniaxial nematic liquid crystals composed of molecules of cylindrical sym- 
metry. A brief summary of the status of theoretical development for the elastic con- 
stants of nematics is presented. Considering a pair potential having both repulsive 
and attractive parts numerical calculations are reported for three systems MBBA, 
PAA and 80CB. For these systems the length-to-width ratio Xq is estimated from 
the experimentally proposed structure of the molecules. The repulsive interaction is 
represented by a repulsion between hard ellipsoids of revolution (HER) and the at- 
tractive potential is represented by the quadrupole and dispersion interactions. From 
the numerical results we observe that in the density range of nematics the contri- 
bution of the quadrupole and dispersion interactions are small as compared to the 
repulsive HER interaction. The inclusion of attractive interaction reduces the values 
of elastic constants ratios. The temperature variation of elastic constants ratios are 
reported and compared with the experimental values. A reasonably good agreement 
between theory and experiment is observed. 

1 Introduction 

In a previous paper [T] (here after reffered to as I), we developed a theory based on 
the density functional formalism j2] for the deformation free-energy of any systems 
with continuous broken symmetry, and applied the theory to derive expressions for 
the elastic constants of nematic and smectic-A phases of uniaxial (-Doo/i) symmetry. 
These expressions of elastic constants are written in terms of order parameters that 
characterize the nature and amount of ordering in the phase and the structural pa- 
rameters involving the generalized spherical harmonic coefficients of the direct pair 
correlation function of an effective isotropic liquid, the density of which is determined 
using a criterion of the weighted density functional formalism |3|. 

In the present paper we restrict ourselves to the uniaxial nematic phase. According 
to the continuum theory ■4J, the bulk elastic properties of nematics can be described 
by three elastic constants, associated with the restoring forces opposing splay (-R'l), 
twist {K2) and bend {K^) distortions. The distortion free energy density is written 
as 

^AAe = \[K^{y ■ Yvf + K2{rv • V X n)2 + ^3(6 x V x n)^] (1) 
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where n, the director, indicates the preferred direction of the long-axes of the molecules. 
These elastic constants play an important role in characterizing liquid crystal displays. 
It is difficult to measure experimentally the absolute values of these elastic constants 
jSlinilZlEl- Their ratios K2/KiSiiad K3/K1 can be measured more accurately [SI 13 EDI- 
For high duty liquid crystal displays [TT], the ratio K3/K1 is desired to be as small as 
possible. The elastic moduli are temperature and density dependent. The dependence 
on the density is pronounced. A number of measurements are reported which show 
that Ki and K2 have weak temperature dependences whereas K3 rapidly increases 
with temperature and that when the nematic-smectic A transition temperature is 
approached from above, Ki does not show any sharp change but K2 and K3 increases 
anomalously. 

The present paper is organized as follows: In sec. 2, we describe, in brief, the 
status of the theoretical development for the elastic constants of the nematics and 
summarize the working equations used in the calculation. The numerical evaluation 
and results are presented in sec. 3 



2 Theoretical development and working equations 

First we shall comment, in brief, on the status of the theoretical development for the 
elastic constants of nematics. A detail of these works are well documented elsewhere 

[1011121 

In the Landau-de-Gennes theory the free-energy density is assumed to be an ana- 
lytic function of the order parameter tensor. To the extent that the order parameter 
is small, the free-energy density is expressed as an expansion in its various orders and 
gradient terms. The elastic free-energy density upto second order in order parameter 
is written in terms of two elastic constants Li and L2 which are related with the 
Frank elastic constants as 

K, = ^-iL^ + h2)-E" (2) 
K2 = ^1,7-2 (3) 

and 

Ki = Ks (4) 
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Thus the elastic constants moduli vary with temperature like P2 . The prediction 
Ki=K3 is not consistent with the experimental observation. However, this is an 
artifact of the derivation in which only the gradient terms of P2 order parameter have 
been considered. If P4 term is also included in the free-energy expansion, all the three 
elastic constants will be different. Using an expansion of the intermolecular potential 
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in spherical harmonics Priest ^3] also arrived at the same conclusion. For the special 
case of dispersion forces, Nehring and Saupe [13] calculated (up to the second order 
P2) the ratio of the elastic constants 

Ki : K2 : K3 = 5 : 11 : 5 

Apart from the temperature dependence of the elastic constants, via P2 order param- 
eter, there is also a dependence on the density. Priest [T3j showed that the deviation 
from the equality Ki = are related in a simple way to the ratio P4/P2 and 

K,/K = A-A'{T\/T^2) + ... (5) 
K2/K = -2A - A'{T\/P2) + ... (6) 
K^/K = A + AA\T\/P'2) + ... (7) 

where K= | {Ki+ K2+ K^). For the case of hard sphero cylinders, the constants 
A and A' were found to depend on the length-width ratio of the molecules. Most of 
the hard rod models ^U] are strictly valid only for the very long and thin rods and 
they usually predict too large a value for K-^ and cannot reproduce the temperature 
dependences of the elastic constants. 

Several workers have evaluated the elastic constants for the van der Walls 
type potential described by hard sphero cylinders with superimposed attractive in- 
teractions. While Stecki and Poniewierski J3] treatment is based on the the direct 
correlation function approach, mean field (MF) approximation has been adopted by 
Kimura et. al. ^|. These authors found that i^s > Ki > K2 and that the tem- 
perature dependence of their ratios is in accordance with the experiment. Using 
generalized van der Walls (GVDW) theory [T7j, which couples the contributions of 
the short-and long-range pair potentials, the elastic moduli were evaluated in a model 
based on the distributed harmonic forces ^H] between the molecules. In this work 
both the repulsive and attractive forces have been considered as distributed along 
the molecules. This approach disregards the temperature dependence of elastic con- 
stants and assumes perfect orientational order i.e. P2 =Pi =1. It is a variant on the 
ideas of Gelbart and co-workers [T7J who studied the combined effect of repulsive and 
attractive forces. Zakhrov [19j evaluated the elastic constants and order parameters 
by using a theory that is based on the method of conditional distribution (MCD) 
pUj . This method introduces a concept of reduced distribution functions which obey 
infinite chains of integro-differential equations. For an arbitrary equations of the 
chains based on the concept of the mean-force potential (MFP) (TT] |221 a truncated 
procedure was adopted. The numerical solution exhibits certain qualitative features: 
The order parameters decrease with increasing volume and temperature. While the 
observed value of Ki/K and K2/K increases with increasing volume, the values of 
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K^/K decreases with it. Ki/K increases strongly with the length- width ratio and 
K2/K decreases with it, and 0.5 < K^/Ki < 3.0 and 0.5 < K2/K1 < 0.8 

Several workers have considered the application of density functional theory to 
study the elastic properties of nematics. A detail account of these works are sum- 
marized elsewhere by one of us ^U]. In this theory exact expression for the elastic 
free- energy is obtained in terms of order parameters and direct correlation functions. 
The correlation of the ordered phases are, in general, not known and hence to be 
approximated. The functional Taylor expansion is performed to obtain the n-particle 
direct correlation function of an inhomogeneous system from the n- and higher-order 
direct correlation functions of an uniform system. The elastic constants of nematics 
were calculated for a number of model systems using approximate forms of the pair 
correlation function of the medium. However, none of these approximate forms of 
the pair correlation function gives the structure of the medium correctly, and so the 
result reported are not expected to be accurate. 

A unified molecular theory was developed by Singh et. al. [1^ for the deformation 
free-energy of ordered molecular phases. This theory is based on the weighted density 
functional formalism [S] and writes the exact expressions for the elastic constants, in 
terms of integrals involving spherical harmonic coefficients of the direct pair correla- 
tion function of an effective isotropic liquid . Adopting the procedure as outlined in 
I, we express the elastic free- energy density for a uniaxial nematic phase of axially 
symmetric molecules. 

J ciri2Q.M(ri2)[(^^)^l^.™(Ax(ri2)) - mi^vu) (8) 

where p„ is the nematic number density and Pi are the Legendre polynomial order pa- 
rameters. The prime on the summation indicates that li and I2 are even fi2 = ri2/|ri2| 
is unit vector along the intermolecular axis, Cg are the Clebsch-Gordan coefficients, 
Ci-^^i^i are the spherical harmonic coefficients of the direct pair correlation function of 
an isotropic liquid, and Ax(ri2)represents the angle between the principal director at 
ri and r2. Confining the variation of director at r2 in a plane, the Yi2m{^x{'''^i2)) is 
expressed in terms of the distortion angle which is assumed to be small. Performing 
the angular integration over fi2 and comparing eq. (8) with eq. (1), the following ex- 
pressions for the elastic constants of uniaxial nematic phase composed of cylindrically 
symmetric molecules are obtained P, 

(3K, = Y.'(3KSi,l2) (9) 
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The explicit expressions for the first few terms of the series can be written as 



(3K,{2, 2) = {^)hplpl[^j,,, - ^J222] (10) 

/3i^i(2,4) = ~{-^)pIP,P,J2,2 (11) 

(3K,{4, 4) = {^)hlPl[V5J,4o - ^>^442] (12) 

PK,{2, 2) = {^)hlPl[lj22o + y|j222] (13) 

/3X2(2,4)=/3Xi(2,4) (14) 

(3K,{4, 4) = (£)^p^p2[^J44o + ^>^442] (15) 

PK,{2,2) ^ PK,{2,2) (16) 

/?i^3(2,4) = -^/?i^i(2,4) (17) 

and 

/3X3(4, 4) = ( A)^p2p2^^_^^^^ ^ _l^^^^^j ^^g^ 

In this expressions the structural parameter Ji-^i^i are defined as 

Jhhi = r^2dri2Ciii^i{ri2) (19) 

3 Calculation and Results 



As obvious from eq. (9) the evaluation of elastic constants requires the values of 
density, order parameters P2 and P4 (as a function of temperature), the structural 
parameters Ji^^i^i ( as a function of density and length-width ratio xq) and potential 
parameters. In I it was shown that the dipole-dipole and dipole-quadrupole interac- 
tions do not contribute to the free energy (and hence elastic constants) of the uniaxial 
mesophases. For the numerical calculation we consider a system of molecules with 
prolate ellipsoidal symmetry interacting via a pair potential, 

U{ri2, ^1, ^2) = {uhER + Uqq + Udis){ri2, ^1, ^2) (20) 

where uher represents the repulsion between hard ellipsoids of revolution (HER) 
and the subscripts qq and dis indicate respectively, the interactions arising due to the 
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quadrupole-quadrupole and dispersion forces. The explicit from of these interactions 
are given in I. 

For evaluating structural parameters , Ji^^i^i, as a function of density and xq, the 
values of spherical harmonic coefficients C/^/j^ for the potential model (20) are needed. 
The C/^/ji can be obtained by solving the Ornstein- Zernike (OZ) equation, using the 
Percus-Yevick (PY) closure relations. As this evaluation is difficult, in a calculation 
only a finite number of spherical harmonic coefficients for any orientation dependent 
function can be handled. It has been shown that for the potential model (20) the 
inclusion of all the harmonics up to indices /i, /2= 4 makes the series fully convergent. 

In this paper we study in detail the influence of HER, qq, dispersion interactions on 
the elastic properties of three uniaxial nematic liquid crystals: p-azoxyanisole (PAA), 
N-p-methoxybenzylidene-p-butylaniline (MBBA) and 4'-pn-octyloxy-4-cyanobiphenyl 
{80CB). For these systems we have estimated the values of length-width ratio xq 
from the proposed structure of molecules as determined by ^'^N nuclear quadrupole 
resonance and line shape studies. The estimated values of xq are 2.68, 2.87 
and 4.88 for MBBA , PAA and 80CB, respectively. The PY closure relation have 
been solved by Ram and Singh j2Hl for the Qher^ huER and Cher harmonics for 
a;o=3.00, 3.25, 3.50 and 4.00. Taking their results we estimated in I C-harmonics 
for the quadrupolor and dispersion interactions. Adopting similar procedure we have 
evaluated here these harmonics for xq = 2.68, 2.87, and 4.88. With known C- harmon- 
ics we evaluated the values of structural parameters as a function of reduced density 
P*( =Pndo) where do is the molecular diameter. In addition to J/j^ai we need the 
values of order parameters P2 and P4 as a function of temperature and of potential 
parameters eo/k, do and quadrupole moment 6. Here eo is constants with unit of 
energy. In the calculation we use the values of P2 as determined by the experiments 
and estimate P4 as P4/P2 — P2 

We calculate the contribution of the individual terms of the series 



Ki = Ki{2,2) + 2K,{2,4)+Ki{4,4) (21) 

for the HER, qq and dispersion interactions for three values of X(){ = 2.68(MBBA), 
2.87(PAA) and 4.88 (80CB)). A number of observations regarding relative contri- 
butions of these interactions terms have been made from this calculation: In the 
density range of nematics the contribution of quadrupole and dispersion interactions 
are small as compared to the repulsive HER interaction. The absolute values of 
elastic constants are sensitive to the values of molecular parameters. The inclusion 
of dispersion and quadrupole interactions reduce the values of elastic constants ra- 
tios. For a given xq, i^rf^^(2,4) is positive where as i^'2^^^(2,4) and K^^^{2,A) 
are negative. Consequently we flnd that K-^^^ > K^^^ > K2^^ and the ratio 
(^K^ER^j^HER-^ > (K^^^/K^^^). This result is in accordance with the simulation 
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work [27]. The absolute values of Kf^^ increase linearly with the temperature and 
the values of their ratios do not change with temperature. In case of the quadrupole 
interaction for a given xq and quadrupole moment B the contribution of -ft'f'^(4, 4)is 
much smaller as compared to Kl'^{2, 2). As the values of increases the contribution 
of each individual terms of the series (21) and Kf'^ increase significantly. The numer- 
ical values of Kl'' and K^'^ are positive where as fTl'^ is negative. The influence of the 
dispersion interaction is small as compared to the quadrupole interaction 

The temperature dependence of the elastic constants are mainly due to the varia- 
tion of P2 and P4 with temperature. We calculate the values of Ki using the experi- 
mental values of P2- Using available experimental data of P2 for PAA |2Zj and MBBA 
PH] we draw a smooth curve and the values of P2 used in the calculation correspond 
to this smooth curve. For 80CB we use the experimental values of P2 as measured 
by Madhusudana and Pratibha [H]. We have found that the values of P4 for PAA as 
estimated from the relation P4 = P2 are in good agreement with experimental data 
P9ij of deuteriated PAA obtained from the coherent neutron scattering experiment. 
So in the calculation we use this estimated values of P4 

Experimentally it is difficult to obtain the absolute values of elastic constants. It is 
the ratios K2/ Ki and K^/Ki which are usually measured more accurately. In addition 
to these ratios the other quantity which one finds accurately from the experiment is 
Ki/ K. In accordance with the experiment we observe that the ratio K^/Ki decreases 
significantly with increasing temperature whereas the ratio K2/K1 is more or less 
independent of the temperature. So in the following figures we show a comparison 
between the theoretical and experimental values of the ratios K^/Ki, Ki/ K and 
K^/K for the PAA, MBBA and 80CB. 

A comparison between the experimental [301 and theoretical values of ratios 
K^l Ki, Ki/ K and K^^/ K for PAA is shown in fig.l. It can be seen that the theoretical 
values are consistent with the experimental data. We have found that the inclusion of 
quadrupole and dispersion interactions decreases the values of elastic constants ratios. 
Figure 1 shows that the ratio K:^/ Ki decreases significantly with temperature. A weak 
temperature dependence is found for the ratios Ki/K. As temperature increases the 
value of K^/K decrease whereas Ki/K increases. As obvious from fig. 2 a similar 
trend in the variation of these ratios has been found for MBBA. In case of 80CB fig. 3 
a similar but more pronounced variation of these ratios with temperature is observed. 
Further, as physically expected near the nematic-smectic A transition a pronounced 
increase in the value of is clearly observed. 

Table 1 shows the variation of elastic constants ratios with length-width ratio 
near nematic-isotropic transition temperature. It can be seen that with xq the ratio 
K^/Ki and K^/K increase whereas Ki/K and K2/K decrease. 

As mentioned in the text, the values of elastic constants are sensitive to the 
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Tabic 1: Ratio of elastic constants of nematic liquid crystals near nematic-isotropic 
transition temperature Tnj 



Xq 


K:i 
K, 


^3 

K 


Ki 
K 


K2 

K 


2.68 


1.211 


1.276 


1.053 


0.747 


2.87 


1.458 


1.34 


0.919 


0.739 


4.88 


1.512 


1.109 


0.917 


0.682 



values of the structural parameters, order parameters and potential parameters. The 
available information about these parameters are not acurate. So as our knowledge 
of these parameters improves, more accurate values of elastis constants will result. 
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Figure 1: Comparison between the calculated (solid lines) and experimental values of 
the elastic constants ratios for PAA as a function of temperature. The experimental 
ratios o[30], *[31] are shown as K^/ Ki ; o and * ; K^/ K ; A and K^/K ; square 
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Figure 2: Comparison between the calculated (solid lines) and experimental values 
[32 ] of the elastic constants ratios for MBBA as a function of temperature. The 
experimental ratios are shown as K^/ Ki : A ; K^/ K : * and Ki/K : o 
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Figure 3: Comparison between the calculated (Solid lines) and experimental [9] values 
of elastic constants ratios for 80CB as a function of temperature. The experimental 
ratios are shown &s K^/ Ki : A ; K^/K: * and Ki/K : o 
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